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I.  INTRODUCTION 

The  T-matrix  or  null  field  method  has  proved  to  be  an  efficient  computational 
scheme  for  analyzing  the  scattering  of  waves  for  several  difficult  geometries.  Its 
broad  applicability  to  many  fields  of  engineering  and  science  in  handling  acoustic, 
quantum  mechanical,  electro-magnetic  or  -static,  and  elastic  problems  has  been 
realized  only  recently,  see  (1-16).  The  T-matrix  approach  has  also  been  extended 
to  treat  increasingly  complicated  situations  such  as  multilayered  scatterers 
(17-25),  multiple  scattering  by  a finite  number  of  scatterers  (19-25),  multiple 
scattering  in  lattices  (26) , and  multiple  scattering  by  statistical  distribution 
of  scatterers  (27).  However,  in  the  theory  of  elasticity  the  analysis  of 
structures  involving  several  layers  with  solid-fluid  interfaces  is  still  lacking. 
The  purpose  of  this  article  is  to  present  a T-matrix  formalism  for  elastic  wave 
scattering  by  multilayered  regions  with  solid-fluid  interfaces.  Here  the  layers 
are  not  necessarily  consecutively  enclosing  and  thus  in  some  regions  there  may  be 
several  different  scatterers .j  The  basic  methods  for  this  extension  are  given  in 
(21-25)  .v>  We  are  using  the  Green's  dyadic  in  an  integral  representation  although 
the  T-matrix  can  be  obtained  without  using  any  Green's  dyadic  (8). 

f 

This  article  consists  of  9 parts.  Here,  in  part  I,  we  introduce  the  basic  problems 
to  be  treated.  In  part  II  we  obtain  the  integral  representations  for  a displace- 
ment field  in  the  case  of  layered  regions.  Then,  we  introduce  the  spherical  wave 
functions  and  use  them  to  expand  the  Green's  dyadic.  In  part  III  we  expand  the 
various  fields  in  spherical  wave  functions.  Using  the  integral  representations 
derived  in  part  II  we  obtain  the  Q-matrix  for  the  outermost  surface  of  a body 
separating  two  elastic  materials.  This  matrix  and  the  T-matrix  associated  with  the 
body  inside  the  next  outermost  surface  form  the  total  T-matrix  for  the  body.  The 
procedure  in  part  IV  is  similar  to  part  III;  however,  we  must  now  deal  with  a fluid 
outside  and  an  elastic  material  inside  the  outermost  surface.  This  new  situation 
requires  the  defining  of  three  matrices  (the  P,  Q,  and  R matrices)  for  the  outer- 
most surface,  in  contrast  to  the  simpler,  "single  surface  matrix"  case  of  part  III. 
These  three  matrices  and  the  T-matrix  associated  with  the  body  inside  the  next 
outermost  surface  form  the  total  T-matrix  for  the  body.  In  part  V we  have  an 
elastic  material  outside  and  a fluid  inside  the  outermost  surface  of  the  body. 

Here  we  have  to  Introduce  the  U,  V,  W and  Y matrices  associated  with  the  outermost 
surface.  These  four  matrices,  together  with  the  T-matrix  associated  with  the  next 
outermost  surface  of  the  two  layered  body,  form  the  total  T-matrix  for  the  body. 


In  part  VI  we  study  scattering  from  several  nonoverlapping  bodies.  In  part  VII  we 
use  the  results  from  the  parts  1I-VI  to  get  a general  procedure  for  obtaining  the 
T-matrlx  associated  with  a multilayered,  finite  body.  These  layers  may  be  either 
elastic  materials  or  fluids,  and  need  not  be  successively  enclosing.  In  part  VIII 
we  review  the  terminology  for  scattering  data.  In  part  IX  we  give  numerical 
results  for  rotatlonally  symmetric  bodies  In  a fluid.  The  wavelengths  under  con- 
sideration are  In  the  longwave-  and  resonance-regions.  We  also  consider  complex 
wave  numbers,  the  moduli  of  which  are  less  than  the  Inverse  power  of  some  charac- 
teristic body  dimension.  The  results  are  presented  In  the  form  of  frequency 
dependent  plots  for  the  blstatlc  differential  scattering  cross  section,  the  total 
scattering  cross  section,  and  the  absorption  cross  section.  We  also  present  some 
polar  plots  of  the  differential  scattering  cross  section. 


II.  INTEGRAL  REPRESENTATIONS  AND  SPHERICAL  WAVE  FUNCTIONS 


Consider  an  Infinite,  homogeneous  elastic  material  with  density  po  and  Lami 
parameters  Xg  and  pg.  A finite  layered  body  Is  Introduced  In  this  material.  The 
body  Is  such  that  between  the  outer  bounding  surface  and  an  Inner  surface  S2 
the  density  Is  p^  and  the  Lami  parameters  X^  and  p^.  The  material  parameters 
Inside  the  surface  S2  can  be  varying  with  the  space  coordinates  In  an  arbitrary 
way.  The  region  outside  the  surface  and  the  region  between  the  surfaces 
and  S2  are  called  region  0 and  region  1,  respectively.  We  assume  that  an  origin 
can  be  chosen  Inside  the  two  surfaces  In  such  a way  that,  for  each  surface,  the 
radius  (from  che  origin)  to  that  surface  Is  a continuous  single-valued  function  of 
the  spherical  angles.  Furthermore,  we  assume  that  both  surfaces  satisfy  the 
Requirement’s  of  the  divergence  theorem.  Let  a monochromatic  displacement  field 
u^,  with  time  dependence  exp(-iwt) , Impinge  on  the  body.  The  scattered  field  ug 
and  the  Incoming  field  u^j  give  the  total  field  in  region  0 as: 


■»  -*1  . -*s 

u0  “ u0  + u0 


(2.1) 


The  total  field  between  the  surface  S]  and 
the  exception  of  source  points  we  have  the 
uA(r) : 


S2  (region  1)  Is  called  uj_.  Then  with 
following  equations  of  motion  for 
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For  the  same  regions  we  define  the  Green's  dyadic  £^(k^,K^, | r-r ' | ) as  the  outgoing 
solution  of 
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W*^  - K1Vx(VxS1)  + - - Ts(r-r')  , (2.3) 

ki 

where  i ■ 0 or  1 and  T is  the  unit  dyadic  and  5 is  the  Dirac  distribution.  The 
stress-dyadlcs  and  trladlcs  for  the  same  regions  are  given  by  (2. A)  and  (2.5), 
respectively. 


(2.4) 
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Ci  " fl,°i 


(2.6) 


The  equations  (2.2)  and  (2.3)  can  be  rewritten  in  terms  of  the  stress-dyadic  and 
triadic,  respectively. 


^’°i^ui^  + pi“  ui  " 0 


(2.2') 


7,£^(£j)  + ■ - T ^ 6 (r-r ' ) 

it, 


(2.3') 


These  latter  forms  may  be  better  suited  for  the  derivation  of  the  integral 
representation.  Here  equation  (2.2')  is  expressed  in  the  same  terms,  which  appear 
in  the  boundary  conditions  of  the  integral  representation.  However,  the  earlier 
forms  ((2.2)  and  (2.3))  may  be  better  suited  for  obtaining  the  basis  functions  or 
Green's  dyadic.  Here  we  can  clearly  see  how  to  separate  the  curl  free  and 
divergence  free  parts.  By  using  equations  (2.2),  (2.3),  and  the  divergence 
theorem  we  can  obtain  the  Integral  representations  for  the  various  regions  (9). 
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Uq(t),  r outside  S^. 


0,  r inside  S^. 


(2.7) 
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r outside  S, 


or 


inside  S2 . 


Here,  the  + and  - subscripts  stand  for  the  limit  values  calculated  from  the 
outside  and  the  inside  of  the  surface,  respectively.  If  the  material  in  one  of 
the  regions  was  a fluid  (u-i  • 0 for  some  i) , we  would  get  (2.9)  and  (2.10)  instead 
of  (2.2)  and  (2.3). 


77*1^  + - 0 


(2.9) 


...  * 


The  only  change  In  the  Integral  representations  (2.7)  and  (2.8)  is  that  we  have  to 
change  the  to  a k^. 

Alternatively,  when  dealing  with  a fluid  we  can  work  with  a scalar  potential 
which  is  related  to  the  displacement  by: 


«i  - • (2.11) 

The  potential  $*(?)  and  the  corresponding  Green's  function  gi(ki, | r-r’ j ) have  the 
following  equations  of  motion: 

(V2+k2)$t  - 0 (2.12) 

(V2**2^  - - <5 (r-r')  . (2.13) 

From  the  equations  (2.12)  and  (2.13)  we  can  derive  the  following  integral  repre- 
sentations valid  only  for  fluids. 


♦£<*>  ♦ ,M  " 

S1 
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4>Q(r),  r outside 
0,  r inside  S^. 


(2.14) 
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inside  S2 . 


(2.15) 


Here  $q,  the  total  field  in  region  0,  is  given  in  terms  of  the  incoming  field  $g 
and  scattered  field  by: 


♦n  “♦£+♦£• 


(2.16) 


To  uae  the  above  Integral  representations  (2.7),  (2.8),  (2.14)  or  (2.15)  one  must 
introduce  boundary  conditions.  These  may  be  expressed  in  terms  of  physical 
observables  such  as  displacement  and  tension  at  the  surface  S^.  The  surface 
tension  ?£  is  expressable  in  a general  material  as  a function  of  the  displacement 
u as  seen  from  equations  (2.4)  and  (2.6).  In  the  particular  case  of  a fluid  the 
surface  tension  can  be  expressed  as  a function  of  the  potential  $q  as: 

C0+  " “ k0X0*0+ftl  P0W  *0+Al  * 


(2.17) 


In  Che  case  of  a fluid,  either  of  che  equations  (2.7)  or  (2.14)  as  well  as  eicher 
of  the  equations  (2.8)  or  (2.15)  can  be  derived  from  the  other. 


It  is  our  aim  to  develop  a matrix  representation  of  the  fields.  The  spherical 
wave  solutions  (2. 18)— (2. 20)  to  eq.  (2.2)  are  well  suited  as  a basis  in  such  a 
representation. 
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Pm(cos8)) 
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i cosm$ 

Y Pm(cos0)}  . 
mn  n | s inm$  , 


(2.19) 


(2.20) 


(2.21) 


where  sq  • 1 and  em  ■ 2 for  m > 0,  e and  o stands  for  even  resp.  odd,  hn  is  the 
(outgoing)  spherical  Hankel  function,  P®  is  the  as^ociat^d  Legendre  polynomial. 

We  are  going  to  use  various  abbreviations  = v*p  = • Further,  the  symbols 

Ou  and  Re  will  represent  the  outgoing  and  regular  functions,  respectively.  For 
example,  OuiJj  ■ ifqt  but,  in  contrast,  Reyq  means  that,  instead  of  using  hn,  we  use 
the  regular  function  jn  (at  the  origin).  This  is  not  necessarily  the  same  as 
taking  the  real  part  of  The  spherical  wave  solution  to  eq.  (2.12)  is: 


>*<r)  " Y™n<0'*)hn<ktr> 
omn  cram  n 1 


(2.22) 


3/2 


When  dealing  with  a fluid  we  let  iji^p  ■ i^p  • 0 and  remove  the  factor  (k^/ic^) 
4>3p . Given  the^above  special  conditions  we  can  interpret  the  expansion  of  the 
Green's  dyadic  G^  in  spherical  wave  functions  for  all  cases  as: 


in 


^i0c1.»c1.|r-r’ |)  - i^(r>)Re^(r<)  , 


(2.23) 


where  r>  and  r<  stands  for  Che  greater  resp.  smaller  of  r and  r'.  The  Green's 
function  g^  can  be  expanded  as: 


g1(k1,|r-T' |)  - ik1^p(r>)Re«»p(r<)  . 


(2.24) 


III.  LAYERED  BODIES  WITH  ELASTIC  MATERIAL  ON  BOTH  SIDES  OF  THE 
BOUNDING  SURFACE 


Consider  an  elastic  body,  with  an  Inclusion,  situated  in  an  infinite  elastic 
material.  Outside  the  bounding  surface  as  well  as  between  and  an  inner 
surface  S5  Che  material  parameters  are  constant  and  given  in  the  notations  from 
section  II.  In  section  II  we  introduced  the  spherical  vector  wave  functions  as  a 
basis  and  expanded  Che  Green's  dyadic  in  this  basis.  Here  we  expand  the 


displacement  fields  In  the  same  basis.  The  expansion  of  the  Incoming  and 
scattered  field  (in  the  Infinite  material)  is  given  by  (3.1)  and  (3.2), 
respectively. 


Here,  r Is  the  length  of  the  radius  vector  from  origin  to  the  closest  singular 
ity  of  the  sources  of  the  incoming  field.  The  symbol  rj;av  is  the  radius  of  the 
sphere  Xi»  with  center  at  the  origin,  circumscribing  Si  (see  Fig.  1).  By  using 
the  integral  representation  (2.7)  and  considering  r outside  the  sphere  X^, 
applying  the  expansions  (3.2)  and  (2.23),  we  get  (3.3),  after  identifying 
coefficients. 


Let  Zi  be  th<j  sphere,  with  center  at  the  origin,  inscribed  in  (see  Fig. 

Considering  r inside  Z,  we  get  (3.4),  in  a similar  way  as  above,  from  (2.7) 


The  total  T-matrix  (32-34)  for  the  layered  body,  defined  by  (3.5),  can  not  yet  be 
obtained  because  we  have  not  yet  specified  the  properties  of  the  body. 


The  relations  (3.3)  and  (3.4)  together  constitute  a formal  definition  of  the 
T-matrix  through  the  two  mapings  from  the  surface  field  to  the  scattered  field  and 
the  incoming  field,  respectively.  To  realize  these  mappings  we  introduce  the 
boundary  conditions  at  Si  and  expand  the  field  on  the  inside  of  Si.  The  welded 
boundary  conditions  (3.6)  and  (3.7)  relate  the  limit  values  of  the  field 
quantities  from  each  side  of  the  surface  S.. 


The  axpanslon  of  the  field  on  the  inside  of  S^  is  given  by 

V " 5(oJRe^J  + 3^q}  * 


The  completeness  of  expansions  of  this  kind  has  been  treated  in  (28-31).  Using 
(3.6),  (3.7),  and  (3.8)  in  (3.3)  and  (3.4),  we  get: 


(ReReQ 


+ ReOuQ 


and 


a ■ - il, (OuReQ1  .o1,  + OuOuQ1  , 01, ) 

q q qq  q qq  q 


(3.10) 


where  Che  matrix  Q is  given  by 

,3 


jOu\JOul-.i+l  Ki  f r-f  /( Ou\-fi\  (OuWi+l 

(Re) (Ref ^qq'  2 J I eillReKrW*q'  “ 

sr+1  u ' ' 

)]' 


(3.11) 


ds  for  i ■ 0. 


Here,  the  first  jRg{  in  the  left  member  of  equation  (3.11)  corresponds  to  the 
selection  of  either  the  regular  or  outgoing  functions  associated  with  the  first 
index  in  Q . Similarly,  the  second  {Rg}  is  related  to  the  second  index.  In 
vector  matrix  notation  we  rewrite  (3.9)  and  (3.10)  as: 

? 


iReReQ1^1  + iReOuQ1!1 
- iOuReQ1^1  - iOuOuQ1]?1 


(3.12) 


(3.13) 


If  the  body  under  consideration  is  homogeneous  (no  inclusion)  which,  in  turn, 
implies  that  r = 0 we  can  solve  the  equations  (3.12)  and  (3.13)  by  elimination 
a^.  In  this  special  case  we  get  the  T-matrix  for  a homogeneous  elastic  body  as: 


of 


- ReReQ1 (OuReQ1)-1 


(3.14) 


However,  if  the  body  has  an  inclusion,  with  bounding  surface  S2»  we  have  to 
Introduce  the  integral^representation  (2.8).  By  using  the  integral  representation 
(2.8)  and  considering  £ outside  the  sphere  we  get  (3.15).  Further,  by  using 
(2.8)  and  considering  r inside  the  sphere  Z2  we  get  (3.16).  The  sphere  Z2  is 
defined  for  the  surface  S2  in  a way  similar  to  the  defining  of  the  sphere  Z^  with 
respect  to  the  surface  (see  Fig.  1). 


Z , (ReReX1  , c.1  + ReOuX1  , B1 , ) 

q qq  q qq  q 

- ReiT1*t1+(u1)]ds  - 0 , 

Z .(OuReX1  a1  + OuOuX1  .B1,) 

q qq  q qq  q 

- ^-V^)  ]ds  - 0 , 


S2 


1+ 


,/‘V^ 


*1+ 


(3.15) 


(3.16) 


where 


(3.17) 


ChtKV  • jp 


'i*  “i 

ds  for  i ■ 1. 

It  Is  a scraighc  forward  matter  to  show  (3.18) 


\Kt 


0 for  OuOu  or  ReRe. 

..  f(+l)  for  0uRe* 

qq'11  1>(-1)  for  ReOu. 


(3.18) 


Expressions  similar  to  (3.18)  has  been  studied  for  the  acoustic  and  electromag- 
netic cases  in  (21)  and  (22),  respectively.  In  (7)  is  expressed  as  the 

difference  between  the  Q-matrices  for  a cavity  and  a rigid  body  and  the  result 

(3.18)  is  obtained  there.  More  interesting  is  perhaps  the  possibility  to  develop 
matrix  formulations  without  using  Green's  functions  or  dyadics.  This  is  done  by 
constructing,  for  each  kind  of  problefn  (acoustic,  electromagnetic  or  elastic),  the 
quantity  corresponding  to  the  integrand  in  (3.17),  see  (8).  Using  the  remark 
about  the  T-matrix,  Just  after  equation  (3.5),  and  the  relation  (3.18)  we  get 

(3.19)  from  (3.15)  and  (3.16). 


I1  - T2^1  (3.19) 

2 

Here  T is  the  T-matrix  for  the  body  inside  the  surface  S2  (i.e.  the  inclusion). 

We  notice  that  the  material  parameters  inside  the  surface  S2  can  be  dependent  on 
the  space  coordinates.  Furthermore,  the  surface  S2  might  be  interpreted  solely  as 
a mathematical  abstraction  (i.e.,  it  doesn't  have  to  be  the  boundary  between 
different  materials).  Except  for  possible  separated  inclusions  within  the 
material  contained  in  and  S2  may  be  identical.  From  (3.12),  (3.13)  and  (3.19^ 
we  derive  the  T-matrix  for  the  entire  layered  elastic  body  as: 


T - - [ReReQ1  + ReOuQ1!2] [OuReQ1  + OuOuQ^T2]-1  . (3.20) 

We  will  end  this  section  by  deriving  the  field  between  the  surfaces  and  S2  in 
the  body.  By  using  the  integral  representation  (2.8)  and  considering  t in  the 
region  between  the  sphere  and  the  sphere  X2,  applying  the  expansions  (3.8)  and 
(2.23),  we  get: 

u - - i Z , (OuReX1  .Re^a1,  + OuOuX1  .Re^1^,)  + 

1 q.q  qq  q q qq  q q 

3 (3.21) 

+ {(*  /*t VRe?J)  ‘“l+  - 

V PjW  b'  ' 

Using  the  remark  about  the  T-matrix,  just  after  equation  (3.5),  and  the  relations 
(3.18)  and  (3.19),  we  get  (3.22)  from  (3.21). 

£3  - Jj(a^Re<£*  + 8^*) 


(3.22) 


This  relation  is  of  course  only  valid  when  the  sphere  X-,  is  truly  inside  the 
sphere  2^. 


IV.  LAYERED  BODIES  WITH  AN  ELASTIC  MATERIAL  INSIDE  THE  BOUNDING 
SURFACE  AND  A FLUID  OUTSIDE 

Consider  an  elastic  body,  with  an  inclusion,  situated  in  an  infinite  fluid. 
Outside  the  bounding  surface  as  well  as  between  S]  and  an  inner  surface  S2  the 
material  parameters  are  constant  and  given  in  the  notations  from  section  II. 
(notice  that  uo  ■ 0)*  The  main  difference  from  section  III  is  that  in  this  sec- 
tion we  only  have  to  deal  with  the  curl  free  basis  functions  outside  (in  the 
fluid).  This,  in  turn,  will  lead  to  some  minor  complications.  The  expansions  of 
the  incoming  and  scattered  field  (in  the  fluid)  is  given  by  (4.1)  and  (4.2), 
respectively. 


Ug(r)  - IapRev^p(r)  for  r < rso 


(4.1) 


u*(r)  - Sf  (r)  for  r > r* 
0 p p 3p 


max 


(4.2) 


so  1 

Here,  r and  rmav  has  the  same  meaning  as  in  section  III.  By  a procedure  similar 
to  the  one  used  in  section  III  we  get  (4.3)  and  (4.4).  Notice  that  we  have  Rq 
here  instead  of  <q  as  in  section  III. 


f - i 
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The  total  T-matrix  for  the  layered  body  is  defined  directly  by  (4.5)  and 
indirectly  by  (4.3)  and  (4.4). 

f - I,T  ,a  , 

P P PP  P 

The  boundary  conditions  at  the  surface  is  given  by  (4.6),  (4.7)  and  (4.8) 


(4.5) 


Al’u0+  " fll*ul-  <for  reSi^ 


(4.6) 


fll*c0+  " (f°t  reSi) 


^tan  " [tl-]tan  * 0 (for  rcSl> 


(4.7) 

(4.8) 


Here,  l )tan  stands  for  the  tangential  component.  The  tangential  components  of 
$04.  and  iS^_  are  unrelated.  However,  because  of  the  relation  (4.9)  we  notice  that 
no  information  about  the  tangential  part  of  or  is  needed. 


Using  (4.6),  (4.7),  (4.9)  and  the  expansion  (3.8)  in  (4.3)  and  (4.4),  we  get: 


fp  " VReReVV  + Re0uQPq,8q,) 


a - - il, (OuReQ  ,a*,  + OuOuQ  .6*.) 

p q pq  q pq  q 


(4.10) 

(4.11) 


where  the  matrix  Q is  given  by 
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OuWl 

fWV 


{2}{£}w  ■ 7?  s/[vfc^-{ 

P0W  ^ 

In  vector  matrix  notations  we  rewrite  (4.10)  and  (4.11)  as: 
t - iReReQa1  + iReOuQ^1 

a " - iOuReQa^  - iOuOuQ^  . 


(4.12) 


(4.13) 


(4.14) 


Notice  that  the  Indices  in  the  matrix  Q belong  to  different  sets  (i.e.,  it  can  not 
be  made  square  by  choosing  different  finite  truncations  in  the  two  indices  n and 
n' . The  most  obvious  counterexample  is  a sphere.)  This  means  that  even  if  the 
body  under  consideration  is  homogeneous  (no  inclusion)  which,  in  turn,  implies 
£hat  r = 0 we  can  not  yet  solve  the  equations  (4.13)  and  (4.14)  by  elimination  of 
a^.  To  overcome  this  we  introduce  the  expansion  (4.15)  of  the  field  on  the  outside 


of  S . 


0+ 


Zy  Reijl? 
PP  3p 


(4.15) 


By  using  the  integral  representation  (2.8)  and  considering  r outside  the  sphere  X,, 
applying  the  expansion  (4.15),  and  applying  the  boundary  conditions  (4.6),  (4.7),^ 
and  (4.8),  we  get  (4.16),  after  identifying  coefficients. 


Z.ReReP  ,y  , + Z, (ReReR  ,a\  + ReOuR  ,6^,) 

p qp  p q qq  q qq  q 


* r 
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(4.16) 


where 


Mv  ■ 77 


(4.17) 


and 


, in  ari  ii  - I-  - 


(4. IS) 
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We  rewrite  (4.16)  in  a vector  matrix  notation  using  (3.19)  and  the  remark  about 
the  T-matrix,  just  after  equation  (3.5). 


ReRePy  + ReReRa1  + ReOuRT2^1  + iT2^1  - 0 
We  also  rewrite  (4.13)  and  (4.14)  after  using  (3.19). 

t - iReReQa1  + iReOuQT2^1 


a * - iOuReQct^  - iOuOuQT“a^ 


(4.19) 


(4.20) 


(4.21) 


If  the  body  under  consideration  is  homogeneous  (no  inclusion)  we  have  T“  = 0.  The 
T-matrix  for  the  homogeneous  body  with  no  inclusion  is  given  by  (4.22). 


T - - ReReQ (ReReR) _1ReReP [OuReQ (ReReR) -1ReReP ] -1 


(4.22) 


The  T-matrix  for  the  entire  layered  elastic  body  with  an  inclusion  possibly  as 
general  as  described  in  section  III,  is  given  by  (4.23). 


T - - (ReReQ  + ReOuQT2) (ReReR  + ReOuRT2  + iT2)_1ReReP < 
x[ (OuReQ  + OuOuQT2) (ReReR  + ReOuRT2  + iT2) -1ReReP ] _1 


(4.23) 


V.  LAYERED  BODIES  WITH  A FLUID  INSIDE  THE  BOUNDING  SURFACE  AND  AN 
ELASTIC  MATERIAL  OUTSIDE 


Consider  a bounded  volume  of  fluid,  with  an  inclusion,  situated  in  an  infinite 
elastic  material.  Outside  the  bounding  surface  S^  as  well  as  between  S]_  and  an 
inner  surface  S2  the  material  parameters  are  constant  and  given  in  the  notations 
from  section  II  (notice  that  - 0) . We  can  proceed  as  in  section  III  using  the 
integral  representation  (2.7),  the  expansions  (2.23),  (3.1),  and  (3.2)  to  obtain 
the  formulas  (3.3)  and  (3.4),  which  are  applicable  here  also.  The  boundary  condi- 
tions are  given  by  (4.6),  (4.7),  and  (4.8)  in  section  IV.  We  now  introduce  the 
expansions  (5.1)  and  (5.2)  for  the  field  on  the  outside  of  and  on  the  inside  of 
S^,  respectively. 


CH- 


Ey  ReiJ^ 
q'q  q 


(5.1) 


- Efa^Rr*1 


e*-  + bJ#  ) 

p 3p  p 3p 


(5.2) 


Notice,  that  outside  (in  the  elastic  material)  we  need  both  curl  and  divergence 
free  basis  functions.  But,  inside  (in  the  fluid)  we  need  only  the  curl  free 
basis  functions.  Using  the  boundary  conditions  (4.6),  (4.7),  and  (4.8)  together 
with  the  expansions  (5.1)  and  (5.2)  in  (3.3)  and  (3.4),  we  get: 


t - iHeReUy  + IReRevJ1  + iReOuVj1 

a - - iOuReUy  - iOuReVa1  - iOuOuvl1  , 
where  the  matrices  U and  V are  given  by: 


(5.3) 


(5.4) 


' ^2  sA({SK)-[R^']tand3 

SKIS'.'  ■ 4 ^o(fcK)-«A-fcR. 

“ {ReK^l^l^lRe^p']'12  * 


(5.5) 


(5.6) 


Notice  that  U is  "square"  but  V is  not.  Here  as  in  section  IV  we  can  see  that 
even  if  the  fluid  is  homogeneous  (no  inclusions)  which,  in  turn,  implies  that 
51  = 0 we  can  not  yet  solve  the  equations  (5.3)  and  (5.4).  We  proceed  in  a way 
similar  to  section  III  with  two  exceptions.  First,  unlike  section  III  where  x^ 
was  used,  we  shall  substitute  k^  in  the  integral  representation  (2.8).  Second, 
unlike  section  III  where  both  curl  and  divergence  free  basis  functions  were  used 
we  shall  expand  the  various  fields  between  and  So  using  curl  free  basis  func- 
tions alone.  By^using  the  integral  representation  (2.8)  (with  k^  instead  of  x^) 
and  considering  r outside  the  sphere  Xi,  applying  the  expansion^  (5.2),  we  get 
(5.7),  after  identifying  coefficients.  Similarly,  considering  r inside  the  sphere 
Z2  we  get  (5.8). 


V(ReReX3p3p'V  + Rfi0uX3p3p»®J’) 
k3  r 

- J[  t^Re^-u^  - Re^p-t^u^lds  - 0 
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Plu  S2 

Notice  that  two  modifications  have  been  made  when  changing  from  an  elastic  medium 
as  in  section  III  to  a fluid.  The  first  is  the  change  from  x^  to  k^  in  the 
integral  representation  (2.8).  The  second  modification  is  the  removal  of  the 
factor  (k^/x^)3'2  in  the  curl  free  basis  function.  These  two  changes  cancel  each 
other  out.  For  this  reason  we  obtain  a subpart  of  the  matrix  X1  as  defined  in 
(3.17).  By  using  the  properties  of  the  matrix  X1  given  by  (3.18)  and  the  remark 
about  the  T-matrix  after  equation  (3.5),  we  get: 


l1  - T2 


(5.9) 


where  Is  Che  T-matrix  for  the  body  inside  the  surface  S2  (i.e.  the  inclusion). 
The  equations  (5.3),  (5.4),  and  (5.9)  are  still  not  sufficient  to  give  us  the 
T-matrix  for  the  total  body.  We  have  to  use  the  integral  representation  (2.8), 
but  now  with  the  expansion  (5.1)  introduced  through  the  boundary  conditions  at  S^. 
By  using  the  integral  representation  (2.8)  (with  kj^  instead  of  k^.)  and  considering 
? outside  the  sphere  X^,  applying  the  expansions  (5.1)  and  (5.2),  and  applying  the 
boundary  conditions  (4.6)  we  get  (5.10). 


E,ReReWwyq,  - g,  (ReReY^.aJ,  + *«Ou*3p3p.eJ,) 
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(5.11) 


(5.12) 


Using  (3.19)  and  the  remark  about  the  T-matrix,  Just  after  equation  (3.5)  we  get 
(5.13)  from  (5.10) . 


ReReWy  - ReReYa1  - ReOuYT2^1  + iT2aL  - 0 
We  also  rewrite  (5.3)  and  (5.4)  after  using  (3.19). 

t - iReReUy  + iReReVa1  + IReOuVT2?1 

a - - iOuReUy  - iOuReVa1  - iOuOuVT2^1 


(5.13) 


(5.14) 


(5.15) 


If  the  body  under  consideration  is  homogeneous  (no  inclusion)  we  have  f H 0.  The 
T-matrix  for  the  homogeneous  Mf luid-body,t  with  no  inclusion  is  given  by  (5.16). 


T - - (ReReU  + ReReV(ReReY)~1ReReW) (OuReU 
+ OuReV (ReReY) -1ReReW) -1 


(5.16) 


The  T-matrix  for  the  entire  layered  "fluid-body",  with  an  inclusion  possibly  as 
general  as  described  in  section  III,  is  given  by  (5.17). 


T - - [ReReU  + (ReReV  + ReOuVT2 ) (ReReY  + ReOuYT2  - iT2)_1ReReW] » 

( 

•[OuReU  + (OuReV  + OuOuVT2) (ReReY  + ReOuYT2  - lT2)"1ReReW]~1 


(5.17) 
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VI.  MULTIPLE  SCATTERING 

Consider  N elastic  bodies  situated  in  an  Infinite  elastic  material.  Outside  the 
bounding  surfaces  of  the  N bodies  the  material  parameters  are  constant  and  given 
in  the  notations  from  section  II.  The  integral  representation  (2.7)  can  be  used 
after  changing  the  region  of  Integration.  We  then  have  the  following  expression: 


ic 

50(?)  + l "^2  • VV^0]ds'  " 


pou  Si 


(6.1) 


uQ(r),  r outside  for  all  i < N. 

0 , r Inside  for  a specific  1 $ N. 

For  the  expansion  of  the  Green's  dyadic  and  the  incoming  field  we  use  (2.23)  and 
(3.1),  respectively.  However,  for  the  scattered  field  we  now  must  deal  with  two 
different  regions.  One  region  is  the  volume  outside  of  the  smallest  sphere 
circumscribing  all  the  surfaces  (see  Fig.  2).  The  other  region  is  the  volume 
inside  the  largest  sphere  inscribed  within  all  the  surfaces  S^.  Both  spheres  have 
center  at  the  origin  0.  The  radii  of  these  spheres  are  called  rout  and  rre8, 
respectively  (see  Fig.  2).  We  now  expand  the  scattered  fields  as  follows: 


«S  - £"  r>c°ut 


(6.2) 


(6.3) 


Before  proceeding  further  we  require  the  translation  properties  for  the  basis 
functions.  These  are  given  by  the  following  formulas: 


E,R(a)  ,tT^,(r)  for  a < r. 

3 qq  q 

t (r+a)  - 

» k-j 

E,o(a)^,Reif^,(r)  for  a > r, 
Re^(rfa)  ■ E,R(a)  ,Rei^,(r)  for  all  a and  r. 

q q qq  q 


(6.4) 


(6.5) 


The  properties  of  the  matrices  R and  a for  the  acoustic  and  electromagnetic 
problems  are  treated  in  (24)  and  (23).  Here  we  simply  need  the  exterior  product 
of  these  two  kinds  of  representations.  The  idea  is  to  express  the  surface 
integrals  in  terms  of  the  radius  vectors  r^  extending  from  the  origins  0^.  These 
coordinate  systems  are  pure  translations  of  each  other  (see  Fig.  2) . By  consider- 
ing r Inside  the  inscribed  sphere  of  surface  S^,  with  center  at  the  origin  0^,  we 
obtain  (6.6)  by  using  the  expansions  (3.1),  (2.23),  (6.4),  and  (6.5)  and  by 
Identifying  the  coefficients  of  the  regular  functions. 


RC(a.)a  - l1  - I o(-S.+a.)TJ  for  i,j  - 1,2,3... N, 
1 U±  1 J 


(6.6) 


with  the  restriction: 


■a 


(6.7) 
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The  superscript  t in  (6.6)  stands  for  the  transpose  of  R (which  in  turn  is  equal 
to  the  Inverse  of  R) . The  double  prime  on  in  (6.7)  indicates  that  is  on  the 
surface  S* . Here  we  have  used  the  following  expressions  for  the  components  of  the 
vectors  a*  and  fl. 
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(6.8) 


(6.9) 


By  considering  r outside  the  sphere,  with  radius  rout,  mentioned  above,  we  obtain 
(6.10)  by  using  the  expansions  (6.3),  (2.23),  (6.5)  and  by  identifying  coeffi- 
cients of  the  outgoing  functions. 


t°at  - 


(6.10) 


Further,  by  considering  r Inside  the  sphere  with  radius  rreg,  mentioned  above,  we 
obtain  (6.11)  by  using  the  expansions  (6.2),  (2.23),  (6.4)  and  by  identifying 
coefficients  of  the  regular  functions. 

res  - {o  at>p  (6.U) 


We  observe  that  the  T -matrix  Ti  for  the  body  bounded  by 
relation  between  f1  and  a*  as  follows: 


the  surface  gives  the 


t1  - T^1  . (6.12) 

We  thus  obtain  the  following  system  of  algebraic  equations  with  unknown 
coefficients  a*,  which  can  be  eliminated: 

RC(a. )a  ■ a*  - I a(-a,+a,)T,aj  for  i.j  - 1,2,3. . .N  , (6.13) 

1 jH  1 3 3 

with  the  restriction:  |-a^+aj|  > r^'  , 

f°Ut  - ntd^a1  . (6.14) 

■?reg  - Zo(a.)T.ai  . (6.15) 

i 1 1 

The  separation  requirement  is  weaker  than  separability  by  planes  or  nonoverlapping 
spheres.  We  note  that  the  result  (6.13)-(6.15)  is  independent  of  whether  the 
infinite  medium  is  a fluid  (with  pg  - 0)  or  a fully  elastic  material  (with  pg  i 0) . 
The  differences  in  the  derivation  when  the  infinite  material  is  a fluid  is  that  kq 
has  to  be  changed  to  kg  and  only  the  curl  free  basis  function  is  needed  (i.e. 
t ■ 3).  Furthermore,  it  is  clear  the  kind  of  material  Inside  is  irrelevant  as 
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Fig.  5 where  Soi»  S02  and  Sj_2  are  separating  different  materials.  This 
configuration  can.  as  shown  in  (21)  and  (22),  be  treated  as  multiple  scattering 
from  two  bodies.  However,  the  origins  of  the  various  coordinate  systems  involved 
muse  be  chosen  so  that  the  geometrical  constraint  in  (6.13)  is  satisfied.  The 
radius  vector  describing  the  surfaces  also  must  be  single  valued. 

We  can  also  study  multiple  scattering  in  lattices,  with  unit  cells  composed  of 
more  than  one  scatterer,  as  in  (26).  In  this  case  we  need  only  plug  in  the  proper 
matrices  pertinent  to  our  elastic-fluid  problem  in  the  formulas,  as  given  in  (26). 


VIII.  CROSS  SECTIONS 


In  this  section  we  begin  with  wave  scattering  in  fluids.  As  mentioned  earlier 

S -fc.  £ 0 + 

this  can  be  described  by  a scalar  potential.  The  expansion  <$q  (r)  - Jf^J(r) 

($p  is  given  by  (2.22))  of  the  scattered  scalar  field  $$  gives  us  the  following 
asymptotic  behavior  of  the  scattered  field: 
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This  leads  us  to  the  following  definition  of  the  scattering  amplitude  A , which  is 
independent  of  the  radius  r: 
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We  are  going  to  define  the  various  cross  sections  (32,33,34)  in  terms  of  incoming 
plane  waves.  For  this  case  the  cross  sections  can  be  interpreted  physically. 

This  physical  interpretation  will  be  treated  later. 

The  incoming  plane  wave  ^ is  given  by: 

$*(r)  - Ce  - ^a^1(a,6)Re<frp(r)  - 4irC£inYp(a,8)Re$p(r)  . (8.3) 

Here,  the  angles  a and  6 are  the  spherical  angles  for  the  incoming  wave  vector  kg 
(see  Fig.  6).  This  wave  vector  is  given  by: 


k'O  ■ kQ(sina  cosS.  sina  sing,  cosa)  . 

We  define  the  differential  scattering  cross  section  do/dfi  by: 


do (q, 8,6,$)  _ Lim  2 
<10(9, $)  = r 


AS(a,B,e,») 


(8.4) 


(8.5) 


where  dfi(9,$)  - sin9d9d$.  Again,  a and  8 are  the  spherical  angles  for  the  incoming 
wave.  Observe  that  both  the  scattered  field  and  the  scattering  amplitude  As  are 
functions  of  these  angles,  as  shown  in  equation  (8.5).  An  incoming  plane  wave, 
with  direction  defined  by  a and  6,  produces  wave  scattering  in  the  direction  given 
by  9 and  #.  The  amount  of  energy  scattered  is  proportional  to  the  differential 
scattering  cross  section  do/dft. 


The  total  scattering  cross  section  o Is  defined  by: 
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Here,  the  coefficients  fp  In  the  expansion  of  the  scattered  field  depend  on  the 
spherical  angles  a and  8,  which,  in  turn,  define  the  direction  of  the  incoming 
plane  wave.  The  total  scattering  cross  section  is  proportional  to  the  total 
amount  of  energy  scattered  in  all  directions. 


The  extinction  (or  total)  cross  section  og  is  defined  by: 

a (o,8)  = Im  ~ As(a,B,a,8)  . 
e Ko  0 


(8.7) 


Here,  Im  stands  for  "the  imaginary  part  of",  and  C is  the  amplitude  of  the  incom- 
ing plane  wave.  The  extinction  cross  section  is  proportional  to  the  total  amount 
of  energy  lost  from  an  incoming  plane  wave  with  direction  given  by  a and  8. 

Finally,  we  define  the  absorption  cross  section  by: 

ff  ( a, 8)  = a (a, 8)  - a(a,6)  (8.8) 

a e 


This  cross  section  is  proportional  to  the  amount  of  energy  that  is  absorbed  by  the 
body  from  the  incoming  plane  wave.  The  angles  a and  0 give  the  direction  of  this 
incoming  wave . 


Instead  of  using  the  potential  to  describe  the  scattering  in  a fluid  we  can  use 


the  displacement  field.  The  expansion  u^(r)  * EfpiJ/^p(r)  of  the  scattered  dis- 
placement field  u5  gives  us  the  following  asymptotic  behavior  of  the  scattered 
field: 


u?(r)  - I(-i)nf  Y (0,$)  — for  large  r.  (8.9) 

v Kq  p P P * 

We  define  the  scattering  amplitude  X,  which  now  is  a vector  in  the  t direction  and 
independent  of  the  radius  r,  by: 

£(©»$)  • ^ re  ^ u^(r)  . (8.10) 

As  before  we  are  only  going  to  define  the  cross  sections  for  incoming  plane  waves. 
For  a fluid  the  incoming  plane  wave  uq  is  given  by: 

uj!j(r)  - kQDeik0  1 - Iapl(a,8)Re$^p(r)  ■ Eap(a,8)Re$jp(r)  . (8.11) 

Here  ajj^  is  the  same  as  in  (8.3).  We  now  define,  for  displacement  fields,  the. 
differential  scattering  cross  section  do/dfl,  the  total  scattering  cross  section  o, 
the  extinction  (or  total)  cross  section  <Je,  and  the  absorption  cross  section  ofl 
by: 


AriMhiM 


dfl(0,$) 


p-*s  1 2 

Lim  2 I V 

r^»  r | -►!  | 2 

I “o' 


2tt  it 


,(a*6)  = / / S dn 

$■0  0*0 


o (a,8)  = In  ^ A(a,8,a,8) 
e K0 


o (a,0)  s a (a, 6)  - a(a,$)  . 

a e 


(8.12) 


(8.13) 


(3.14) 


(8.15) 


Here  D Is  Che  amplitude  of  the  incoming  plane  wave.  In  section  II  we  defined  the 
relation  uo  ■ V#q.  From  this  relation  and  by  comparison  of  equations  (8.3)  and 

(8.11)  we  get  the  relation  D - ikgC  between  the  amplitudes  for  the  displacement 
field  and  scalar  field.  We  notice  that  the  same  T-matrix  relates  the  incoming 
field  to  the  scattered  field  in  Che  potential  and  displacement  field  descriptions. 
This  fact,  together  with  equation  (8.11)  and  the  relation  D ■ lkQC,  leads  us  to 
the  following  two  relations:  fp  - k<)fp  and  A - ik()As.  From  the  above  we  can  see 
that  the  cross  sections  defined  by  (S.§)-(8.8)  are  the  same  as  those  defined  bv 

(8. 12) -(8. 15) . 


Ue  now  Cum  to  the  case  where  the  unbounded  medium  is  an  elastic  material.  In 
this  case  we  have  to  deal  with  both  the  divergence  free  and  curl  free  basis 

functions,  (i.e.  t ■ 1,  2,  and  3).  The  expansion  u*(r)  ” If  il^(r)  (where 

q = romn)  of  the  scattered  displacement  field  uq  gives  us  the  following  asymptotic 
behavior  of  the  scattered  field: 

\ tl  r*  / \ 

55<?>  * { ^-[(-  + £2ps2p<e'*>)  v- 


(8.16) 


f3pVM) 


LkQr  , 

~ J 


for  large  r. 


The  functions,  ATp  for  x ■ 1,  2 are  purely  transversal.  The  function  ^3p  is 
purely  longitudinal.  In  Che  case,  with  an  unbounded  elastic  medium  we  have  two 
kinds  of  incoming  displacement  fields.  The  first  kind  is  the  purely  longitudinal 
wave  given  by  (8.11)  with  a slight  modification.  Here,  in  the  case  with  the 
elastic  medium,  we  must  introduce  a factor  (kq/^q)^^  before  the  expansion  coeffi- 
cient ap  because  of  a change  in  the  normalization  of  the  curl  free  basis  function 
\Wp  when  shifting  from  a fluid  to  an  elastic  medium.  The  second  kind  of  incoming 
displacement  field  is  the  purely  transversal  one  given  by: 


IJ(r)  - E.lk°'r 


Z a 

x-1,2  TP 

P 


(a,8)ii^p(r) . 


(8.17) 


Here  £ is  the  amplitude  vector,  which  is  orthogonal  to  t.  The  angles  a and  8 are 
the  spherical  angles  for  the  incoming  wave  vector.  We  define,  for  the  displace- 
ment fields,  the  differential  and  total  scattering  cross  sections  by  (3.12)  and 

(8.13),  respectively.  Notice  that  we  can  separate  the  longitudinal  and  the 


transversal  parts  of  the  different  cross  sections, 
section  can,  by  using  the  relation  f X ,ds2  ■ 5 

* Q d 


a(a,6) 


q q qq 

Firtr 2 - 'V2*; 

lu0l  <0 


The  total  scattering  cross 
i,  he  expressed  as: 


(.8.18) 


IX.  NUMERICAL  RESULTS 

Elastic  wave  scattering  (single  scattering)  by  homogeneous  bodies  of  elastic 
materials  or  cavities  on  an  infinite  elastic  material  is  studied  numerically  in 
(12,14).  Multiple  scattering  of  elastic  waves  by  layered  bodies  (several  bodies) 
in  an  infinite  elastic  material  is  studied  numerically  in  (19).  Numerical  results 
for  the  scattering  of  waves  by  a homogeneous  body  in  a fluid  are  given  in  (15). 

For  more  wave  number  dependent  data  see  (16).  In  this  article  we  shall  give  some 
numerical  results  for  homogeneous  and  layered  elastic  bodies  in  a fluid.  The 
surfaces  of  the  layers  in  the  bodies  are  prolate  spheroids  with  axis  ratios 
ai/bi  - 2.0  (see  Fig.  7).  The  ratios  ao/a^  between  the  semiaxes  of  the  inner  and 
outer  surfaces  are:  0.0,  0.25,  and  0.9.  The  inner  surface  always  encloses  a 
cavity.  The  material  parameters  of  the  fluid  and  of  the  elastic  solid  arc  given 
in  Table  1. 

We  choose  the  2-axis  as  the  axis  of  rotational  symmetry.  Without  any  loss  of 
generality  we  set  the  spherical  angle  6 for  the  incoming  wave  (see  Fig.  6)  equal 
to  zero.  This  choice  reduces  the  number  of  matrix  elements  needed  in  the  finite 
truncation  by  a factor  of  0.5.  In  the  diagrams  studied  here  we  shall  consider 
only  incoming  waves  with  the  spherical  angle  a - 0.0.  In  the  polar  plots  we  shall 
study  the  differential  scattering  cross  section  do/dO  as  a function  of  the 
scattered  angle  0.  In  this  context  the  other  scattered  angle  $ will  take  on  the 
values  0.0  or  180.0.  In  the  Cartesian  plots  we  shall  study  the  differential  and 
total  scattering  cross  sections  as  a function  of  kga^.  We  only  consider  the 
scattering  angles  0-90  and  $ - 0 in  these  plots.  For  both  the  polar  and 
Cartesian  diagrams  we  shall  consider  the  influence  of  losses  in  the  elastic 
material.  The  complex  wave  numbers  associated  with  these  losses  are  given  in 
Table  1. 

Gauss-Legendre  quadrature  formulas  were  used  to  generate  the  matrix  elements  of  Q, 
P,  and  R defined  in  Eqs.  (4.12),  (4.17),  and  (4.18).  In  the  homogeneous  case  the 
matrix  R was  inverted  by  Gaussian  elimination.  However,  the  next  inversion  in 
Eq.  (4.22)  was  done  by  Schmidt  orthogonalization  for  the  case  of  loss  less 
materials.  Here  we  used  a computer  routine  developed  by  P.  C.  Waterman  (35) 
wherein  the  symmetric  and  unitarity  properties  of  the  S-matrix  (S  - 1-2T)  have 
been  used  to  optimize  the  inversion  procedure.  In  the  layered,  loss  less  case  we 
could  use  only  the  Schmidt  orthogonalization  process  for  the  T-matrlx  of  the 
cavity.  The  other  inversions  in  this  case  were  done  by  using  Gaussian  elimination. 
In  the  cases  dealing  with  losses  the  S-matrix  is  not  unitary,  and  we  have  to  use 
Gaussian  elimination. 

# 

The  T-matrix  for  rotatlonally  symmetric  bodies  is  diagonal  in  the  azimuth  index. 

The  two  submatrices  (of  the  T-matrix)  which  corresponds  to  the  two  azimuth  indices 
0 and  1 were  given  the  same  dimension  in  the  finite  truncation.  However,  for  each 
unit  increase  in  the  value  of  the  azimuth  index  beyond  1 we  decreased  the  dimen- 
sion of  the  associated  submatrix  by  one.  The  dimension  of  the  two  biggest  sub- 
matrices was  progressively  increased  from  8 at  kga^  - 0.1  to  12  for  k$a^  - 3.2. 


Figures  8-18  are  plots  of  the  cross  sections,  all  normalized  with  respect  to  a|, 
for  a fictitious  solid.  The  properties  of  this  solid  are  given  in  Table  1. 

Figures  8-13  are  polar  plots.  Figures  14-18  give  the  cross  sections  as 
functions  of  kga^.  Figure  captions  give  details  of  the  individual  plots.  The 
spectra  indicate  chat  the  very  sharp  resonances  for  a loss  less  solid  almost 
completely  disappear  with  the  introduction  of  frequency  dependent  complex  wave 
numbers  in  the  solid.  To  obtain  these  sharp  peaks  we  used  a step  size  of  0.025  in 

*c0al* 

The  method  used  in  this  paper  is  to  our  knowledge  the  only  one  working  for  the 
problems  treated  here.  This  method  shares  with  many  other  methods  (in  other 
areas)  the  lack  of  a rigorous  convergence  proof.  The  parameters  we  have  to  play 
with  in  order  to  obtain  numerical  convergence  are;  the  number  of  integration 
points  and  the  dimension  of  the  matrices  in  the  finite  truncation.  We  have  made 
computations  using  the  Gauss  elimination  and  numerically  tested  the  relations 
(9.1)  and  (9.2)  for  a loss  less  solid. 

r - T (9.1) 

I'T  - - Real  part  of  T (9.2) 

Here,  t means  the  Hermite  conjugate.  The  discrepancy  from  these  relations  was 
measured  relative  to  the  maximum  element  of  Che  T-matrlx.  We  also  use  the  ratio 
between  the  absorption  cross  section  and  the  total  scattering  cross  section  as  a 
convergence  measure  for  the  loss  less  solid.  For  homogeneous  bodies  our  measures 
were  usually  better  than  10~3.  However,  for  the  layered  bodies  they  were  about 
2*10~2.  In  one  exceptional  case,  near  a peak,  we  got  as  bad  a measure  as  0.2.  It 
is  clear  that  some  matrices  will  become  more  illconditioned  for  wave  numbers  close 
to  peaks.  It  is  also  clear  that  matrices  with  the  combination  OuOu  are  especially 
difficult  to  obtain  with  good  numerical  accuracy. 
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Fig.  8. a.  Loss  less  solid. 

Fig.  8.  Polar  plots  of  the  differential 
spheroids  with  kga^  ■ 1.0. 


Fig.  9. a.  Loss  less  solid. 

Fig.  9.  Polar  plots  of  the  differential 
spheroids  with  kga^  » 3.0. 


Fig.  10.  Polar  plot  of  differential 
scattering  cross  section 
for  a layered  spheroid  with 
a2^al  “ 0-25  and  kQa^  * 1.0. 


Fig.  12.  Polar  plot  of  differential 

scattering  cross  section  for 
a layered  spheroid  with 
a2^al  * 0.90  and  kQai  • 1.0. 


Fig.  8.b.  Lossy  solid. 

scattering  cross  sections  for  homogeneous 


; Fig.’  9.b.  Lossy  solid. 

scattering  cross  sections  for  homogeneous 


Fig.  11.  Similar  to  Fig.  10  but  here 


Fig.  13.  Similar  to  Fig.  12  but  here 
*0*1  - 3.0. 


Differential  scattering  cross 
section  for  spheroids  for 
a ■ 0,  8 ■ 0,  0 ■ 90  and  <J>  ■* 

0;  for  a2/a^  • 0.0  loss 

less  solid;  - - - for  a->/a^  - 
0.0  lossy  solid;  for 

a2^al  " 0*25  loss  less  solid. 


k0°l 


Fig.  15.  Similar  to  Fig.  14  but  here 
only  for  m 0*90  loss 

less  solid. 


Fig.  16.  Total  scattering  cross  01 

section  for  spheroids  for 
a ■ 0,  8 ■ 0,  0 ■ 90  and 

^ ■ 0;  for  a2/ai  ■ 

0.0  loss  less  solid;  

for  a2/a^  ■ 0.0  lossy 
solid;  for  a2/a^  ■ 

0.25  loss  less  solid.  0 


*0  al 
Fig.  16. 
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